Differential rotation of main sequence F stars 
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Abstract. The differential rotation of a 1.2 Mq zero age MS star (spectral type F8) is computed and the results are compared 
with those from a similar model of the Sun. The rotation pattern is determined by solving the Reynolds equation including the 
convective energy transport. The latter is anisotropic due to the Coriolis force causing a horizontal temperature gradient of 7 
' K between the poles and the equator. Comparison of the transport mechanisms of angular momentum (the eddy viscosity, the 

' L-effect and the meridional flow) shows that for the F star the L-effect is the most powerful transporter for rotation periods of 7 

. d or less. In the limit of very fast rotation the L-effect is balanced by the meridional flow alone and the rotation is nearly rigid. 

The rotation pattern found for the F star is very similar to the solar rotation law, but the horizontal shear is about twice the solar 
value. As a function of the rotation period, the total equator-pole difference of the angular velocity has a (slight) maximum 
• at a period of 7 d and (slowly) vanishes in both the limiting cases of very fast and very slow rotation. A comparison of the 

^0 ' solar models with those for the F-type star shows a much stronger dependence of the differential surface rotation on the stellar 

luminosity rather than on the rotation rate (see Fig.|3}. 
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1. Introduction 

The differential rotation of the solar photosphere has been ob- 
served for a long time. Depending on the method used, the ro- 
tation period at the equator is found to be 20-30 percent shorter 
than at the poles. Over the last 15 years helioseismology has re- 
vealed that this phenomenon is not restricted to the atmosphere, 
granulation, and supergranulation layers but persists through- 
out the whole convection zone. Moreover, the gas rotation does 
not show the expected increase with depth. Instead, the pat- 
tern observed at the surface is found throughout the whole 
convection zone while the core appears to rotate rigidly at the 
same rate as the convection zone does at intermediate latitudes 
(Thompson et al. 2003). 

The solar rotation pattern can be explained with angular 
momentum transport by the gas flow in the convection zone. 
Assuming that Reynolds stress is the only relevant transporter 
of angular momentum, Kiiker et al. (1993) found a remarkably 
good agreement between the results from their model based 
on the theory of the Reynolds stress and the observed rotation. 
The model, however, does not explain the impact of other trans- 
porters, especially the meridional flow. Kitchatinov & Riidiger 
(1995) and Kiiker & Stix (2001) presented a refined model 
that included both the meridional flow and the convective heat 
transport. They showed that the flows driven by the differen- 
tial rotation and a small horizontal temperature gradient due to 
anisotropic convective heat transport essentially have opposite 
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directions, and roughly balance each other in the solar convec- 
tion zone. 

Hall ( 1 99 1 ) found differential rotation for a number of mag- 
netically active stars from the variation of spot rotation peri- 
ods over the stellar activity cycle. The variation of stellar ro- 
tation periods over the activity cycle was also used to measure 
the surface differential of active lower main-sequence stars by 
Donahue et al. (1996). Messina & Guinan (2003) derived sur- 
face differential rotation from photometric data from a moni- 
toring program that studies a sample of stars that resemble the 
Sun in an earlier state of its evolution. 

The number of single main-sequence stars to which 
Doppler imaging has been successfully applied is still small 
(Strassmeier 2002). For sufficiently fast rotators, surface differ- 
ential rotation can be detected, however, from the broadening 
of spectral lines. Reiners & Schmitt (2003a,b) carried out mea- 
surements for F stars with moderate and short rotation periods. 
They found differential rotation to be much more common for 
stars with moderate rotation rates than for very rapid rotators. 

Here we present a model for the differential rotation of a 
main sequence star of spectral type F8. This star, with 1.2 solar 
masses, lies just below the maximum mass for main-sequence 
stars with outer convection zones and therefore represents the 
upper end of the lower main-sequence in the context of differ- 
ential rotation and stellar activity. With a convection zone depth 
of 160,000 km and a surface gravity roughly equal to the solar 
value the main difference from the Sun is the luminosity, which 
is 1 .7 times the solar value. 
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Ihe radiative heat ilux reads 



The basic setup is the same as in Kiiker & Stix (2001) and very 
similar to Kitchatinov & Riidiger (1999). The stellar parame- 
ters are taken from a model of a 1.2 Mq ZAMS star computed 
by Granzer (2002). The fractional depth of the convection zone 
is 20 % of the stellar radius of 1 . 1 3 solar radii. While the depth 
of the convection zone is similar to that of the Sun, the greater 
luminosity enforces a larger convective heat flux. The density 
stratification is steeper than in the solar convection zone, lead- 
ing to smaller values of the density and pressure scale heights, 
and hence smaller values of the mixing-length. As a conse- 
quence, the convection velocity must be larger than in the solar 
convection zone. 

We treat rotation and large-scale meridional flow as a per- 
turbation with the stratification from a standard stellar evolu- 
tion code as the unperturbed state. For rotation periods longer 
than one day the deviation from spherical symmetry is small. 
The unperturbed state is therefore assumed spherically sym- 
metric. 

The equations to be solved are the conservation laws of an- 
gular and linear momentum, and the equation of convective 
heat transport. With the mean-field ansatz, m = m + u' , the 
equation of motion for the mean velocity field u reads 



du 
'dt 



-W-pQ~VP + pg, 
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with the one-point correlation tensor Qij = {u[u'j). We define 
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and split it into its horizontal average and the deviation from 
that average: 



f3{r,9)^l3\r)+f3\r,9), 
where 

f3\r) = - I I3sin0d0. 
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The convective heat diffusion coefficient is defined in terms of 
(3^, using the form derived by Kitchatinov & Riidiger (1999) 
instead of standard mixing-length theory: 



A similar expression holds for the turbulence viscosity, 
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The heat transport is then described by the transport equa- 



tion 



pT^ = -V • (F^™^ + F'-^^) + pcpu ■ (3\ (7) 



with the convective heat flux. 
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with the opacity k. The advection term in Eq. contains 0^ 
instead of /3 to remove the convective instability from the sys- 
tem. The radiative and convective transport terms account for 
the entire radial heat transport. Heat advection by the large- 
scale meridional flow is included only to allow for a back- 
reaction against the baroclinic term, ^(Vp x Vp)^ in Eq. Q. 
Assuming hydrostatic equilibrium, 

VP = -fifp, (10) 



we can rewrite the baroclinic term in terms of the specific en 
tropy: 

9 ds 
rcp do ' 



\(Vp X VF)^ = — ^(Vs X VP)4 
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The tensor Xij contains the anisotropy of the convective heat 
transport (Kitchatinov et al. 1994; Kapyla, Korpi & Tuominen 
2004, Rudiger et al. 2004). 

3. Results 

The model has first been applied to the Sun by Kiiker & Stix 
(2001). With the original value of 4/15 for and a value of 5/3 
for the mixing-length parameter, a value of 0.2 was obtained for 
the normalised horizontal shear, 
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where il^q and fipoic are the rotation rates at the equator and 
the poles. As the resulting horizontal shear is smaller than the 
observed 30 %, we treat the viscosity parameter Ci, as a free pa- 
rameter and vary it to improve the reproduction of the observed 
rotation law. The turbulent heat transport coefficient c^^ then 
results from the requirement that the Prandtl number keeps its 
value of 0.8. Test computations have shown that for the Sun the 
latitudinal shear is 20 % for = 4/15 and 30 % for = 0.15. 
We therefore carry out the computations for our model F star 
with this value. 

3. 1. The rotation law 

Figure \l\ shows the rotation pattern for the F star model with 
Ci, = 0.15 and q;mlt = 5/3 and rotation periods of 4 d, 7 d, 14 
d, and 27 d. We define the mean angular velocity as the ratio of 
the angular momentum and the moment of inertia, i.e. 



/ pr^ sin^ endrde 



J pr^ sin^ 9drd9 ' 
The mean rotation period is then 

P 



(13) 
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which (naturally) corresponds to the rotation period at a lati- 
tude of 30° for a solar-type rotation law. 
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Fig. 1. Normalised rotation rate as a function of the fractional stellar radius for the 0°, 15, 30, 45, 60, 75, and 90° latitude (from 
top to bottom in each diagram) for rotation periods of 4 d, 7 d, 14 d, and 27 d (from left to right). 



Fig. 2. Isocontours of the rotation rate for average rotation periods of 4 d, 8 d, 16 d, and 27 d (from left to right). 



In Fig.n]the rotation rate is plotted vs. radius for latitudes 
at 0° (equator), 15, 30, 45, 60, 75, and 90° (poles). In all cases 
the equator rotates faster than the poles (solar-type differen- 
tial rotation), but the amplitude of the relative shear varies with 
the rotation rate. The model with the fastest rotation yields the 
most rigid rotation law. Note also that for fast rotation all lines 
are practically horizontal, i.e. there is no radial shear, while 
the model with a rotation period of 27 d shows a substantial 
decrease of the rotation rate with increasing radius at high lati- 
tudes. This change of the rotation pattern can also be seen in the 
contour plots shown in Fig.|2] While the isocontours are mainly 
radial for fast rotation the slow rotation case (27 d) shows a 
disk-shaped pattern at high latitudes. Figure |3] shows the total 
horizontal shear, Sil, as a function of the rotation period for the 
F8 star and the Sun for both values of the viscosity parameter 
Cy. All four curves show a decrease of the shear with increas- 
ing rotation period for slow rotation. For rapid rotation, the be- 
haviour is the opposite: The rotation becomes more rigid as the 
period decreases, i.e. the rotation gets faster Between the two 
limiting cases the total horizontal shear assumes its maximum 
value. The period at which the maximum is reached is about 
one month in case of the Sun, and 19 d for the F star. None of 
the curves shows a sharp peak. There is always a broad range 
of nearly constant shear around the maximum. For both stars, 
the shear strongly depends on the viscosity parameter. 



15 




Ol , ... I , , , , ... I 

10 100 
ROTATION PERIOD [day] 



Fig. 3. The total horizontal shear for the F8 star as a function 
of the rotation period. Solid line: Ci, = 0.15, dot-dashed line: 
= 4/15. Dashed line: solar-type star, Ci, = 0.15. Dotted 
line: solar- type star, = 4/15. 



3.2. The meridional flow 

Figure 13 shows the meridional flow patterns. For P = 4 d the 
figure shows only one flow cell per hemisphere, with the flow 
directed towards the poles in the upper part of the convection 
zone and towards the equator at the bottom. For a rotation pe- 
riod of 7 d, a second flow cell is visible. It is directed towards 



Fig. 4. Stream lines of the meridional flow for rotation periods of 4 d, 7 d, 14 d, and 27 d (from left to right). 



the equator both at the top and the bottom of the convection 
zone, but towards the poles in the bulk of the convection zone. 
At a rotation rate of 14 d the upper cell has become larger, and 
for P = 27 d it dominates the flow pattern, with the counter- 
clockwise flow restricted to a thin layer at the bottom. 

In Fig.|5lthe maximum flow speed at the bottom of the con- 
vection zone is shown. The sign refers to the direction of the 
flow at the point where it is fastest. A positive sign means that 
the flow is toward the equator, negative values indicate pole- 
ward flow. For the F8 star and — 0.15 the value of the speed 
decreases from 10.4 m/s for P^ot = 4 d to -6.3 m/s at P = 60 
d. The change of the flow direction at the bottom of the con- 
vection zone should have dramatic consequences for the stellar 
dynamo, if indeed the form of the butterfly diagram is domi- 
nated by the meridional flow at the bottom of the convection 
zone. 

For rotation periods longer than Prot = 60 d the speed in- 
creases as the absolute value decreases. For fast rotation the 
flow is faster in the Cu = 0.15 case than in the = 4/15 case, 
while for slow rotation the latter gives the faster flow. 

3.3. Differential temperature 

Figure |6l shows the variation of the temperature with latitude at 
the upper boundary for the solar-type star and the F8 star. The 
horizontal average has been subtracted. In both cases the polar 
region is slightly hotter than the equator For the G2 star the 
difference is 1.7 K, for its F8 counterpart 7.0 K. The average 
temperatures are 2.8 x 10^ K and 1.57 x 10^ K, respectively. 

4. Discussion 

For the same set of parameters, the F star shows a larger dif- 
ference between the equatorial and polar rotation rates. As the 
convection zones are of roughly the same depth, this suggests 
that differential rotation is determined by the stellar luminosity 
rather than the rotation rate. Both mixing length and convec- 
tive turnover time in the solar convection zone are larger than 
the corresponding values in the F star convection zone. The op- 
posite holds for the convection velocity. The mixing length in 
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Fig. 5. Meridional flow speed in m/s at the bottom of the con- 
vection zone as a function of the rotation period. Positive value 
correspond to gas motion toward the equator Solid line:F8 star, 
— 0.15, dashed line: F8 star, c^, = 4/15. Dashed line: solar- 
type star, Ci, = 0.15. Dotted line: solar-type star, — 4/15. 
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Fig. 6. The temperature variation (in K) at the upper boundary. 
Left: solar-type star, Prot 27 d, right: F8 star, P = 10 d. The 
horizontal average has been subtracted. 



the F star is only slightly smaller than in the solar convection 
zone, by about as much as the convection zone is shallower. 
The convection velocity is twice that of the Sun because of the 
higher luminosity. 
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Fig. 7. The global convective turnover time as a function of 
the rotation periods for the F star and the solar-type star Solid 
line: F star, — 0.15. Dash-dotted line: F star, = 4/15. 
Dashed line: solar-type star, Ci, = 0.15. Dotted line: solar-type 
star, Ciy = 4/15. 

A comparison between the two convection zones can be 
summarised as follows. Since the stellar radii and depths of 
the convection zone are very similar, the mixing lengths are 
roughly the same, too. The higher luminosity of the F star re- 
quires a larger convective heat flux, therefore larger velocities, 
and hence shorter convective turnover times. As a result, the 
difference in the Coriolis number between an F8 star with a 7 
d rotation period and a G2 star rotating with a period of 27 d is 
less than a factor of two, half the ratio of the rotation periods. 

Kim & Demarque (1996) define the global turnover time, 

= / (15) 

as an integral over the entire extent of the convection zone. 
With a value of 1.86 for the mixing length parameter they find 
Tc = 13.95 d for a 1.2 Mq star. Figure0shows this quantity 
as resulting from our model for the F8 star as a function of the 
rotation period. While not constant, Tg is a slowly varying func- 
tion of the rotation period. Rapid rotation enforces shorter time 
scales because the efficiency of the convective heat transport is 
reduced. As both the convective heat flux and the mixing length 
are prescribed by the boundary conditions and the stellar strat- 
ification, respectively, the only way for the system to adjust is 
by increasing the convection velocity, and hence shortening the 
turnover time. The values found are smaller than the 13.95 d of 
Kim & Demarque (1996), but the curve in Fig.Qshows an in- 
crease towards slower rotation. Using standard mixing-length 
theory, their model corresponds to ours in the limit of slow ro- 
tation. Note that the global convective turnover time is about 
twice as long as its locally defined counterpart. 

Our model predicts the surface rotation pattern of lower 
main-sequence stars to depend on luminosity much more than 
on the rotation rate. This finding is in agreement with the 
observed surface differential rotation of AB Dor (Donati & 
Coflier Cameron 1997), PZ Tel (Barnes et al. 2000), and LQ 
Lup (Donati et al. 2000). These stars, though rotating much 
faster than the sun, show surface differential rotation similar to 
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Fig. 8. The lapping time as a function of the rotation period for 
the F star (solid line) and the solar-type star (dashed). 

that of the Sun. LQ Lup (RXJ1508.6-4423) is a post T Tauri 
star of spectral type G2. With a rotation period as short as 0.31 
d the surface shear Sfl = 0.13 d^^ is only about twice the solar 
Sn = 0.07 d^i. AB Dor and PZ Tel ai-e PMS stars of spectral 
type KO and K9/G0, respectively, with rotation periods of 0.5 
d and 1 d. With Sn = 0.056 d^^ (AB Dor ) and Sil = 0.075 
d^^ (PZ Tel) the surface shear values of these stars lie close to 
the solar value, with the more luminous AB Dor also showing 
more surface shear 

Reiners and Schmitt (2003 a,b) found values between 10 
and 30 for the lapping time, Piap = 2tt/ ((557). Figure |8] shows 
the lapping time vs. the rotation period for the F and solar-type 
stars for c^y = 0.15 from our model. For both stars there is little 
variation except for very long periods, where the lapping time 
increases profoundly. The value for the solar-type star is about 
lOOd at the period of maximum shear That of the F star 50 d, 
2-3 times larger than those found by Reiners & Schmitt for the 
same interval of rotation periods. 

Reiners & Schmitt found differential rotation to be much 
more common for F-type stars with moderate rotation rates 
than for the most rapidly rotating ones. As the method they 
use detects differential rotation only when the relative shear, 
Sil/il exceeds certain threshold value, their finding is not sur- 
prising. As the total shear is a slowly varying function of the 
rotation rate, the relative shear is roughly proportional to the 
rotation period, i.e. 6il/il oc Piot The non-detection of differ- 
ential rotation therefore does not rule out its presence for rapid 
rotators and it is actually predicted by our model. 

Several observers have studied the connection between ro- 
tation period and surface rotational shear for late-type stars. 
The result is usually expressed as a power law, i.e. 

sn cx PZt (16) 

In this representation a constant shear reads n 0. Figure |3] 
shows the surface horizontal shear as a function of the rotation 
period for the F8 star and for the solar model. The functions 
Sil — Sil{Prot) are not power laws, especially not for the F8 
star. A value of 0.5 for the exponent n gives a good fit for the 
latter with = 0.15 in the interval from 4 d to 7 d while a 
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Fig. 9. Horizontal flux of angular momentum in the convection 
zone of the solar model (Pjot = 27d, c^, = 0.15) with the baro- 
clinic term included. Solid: contribution of the meridional flow. 
Dashed: L-effect. Dotted: viscosity. Dash-dotted: total angular 
momentum flux. 

smaller value of -0.5 is found for periods between 14 d and 20 
d. Around Piot = 10 d we find n = 0, which corresponds to 
the maximum of 6Vl = (5i7(Pi.ot)- 

Hall (1991) found n = -0.21 while Donahue et al. (1996) 
derived n = —0.7 from observations of the Ca II chromo- 
spheric lines. Schmitt & Reiners report a value of —0.53, and 
Messina & Guinan derive an exponent of —0.6 (see Riidiger 
& Hollerbach 2004). The observed values tend to be smaller 
than the value of 0...0.5 which our model predicts for rapid ro- 
tation. A possible explanation is scatter caused by the variety 
of spectral types in the observed samples. 

The anisotropy of the convective heat transport plays a cru- 
cial part in maintaining the latitudinal shear. The baroclinic 
term that arises from the horizontal temperature gradient coun- 
teracts the shear along the z-axis as a driver of meridional flow. 

Brun and Toomre (2002) have compared the contributions 
of the L-effect, viscosity, and meridional flow to the total hori- 
zontal flux of angular momentum, 

/■top 

/(6») = / Fg{r,d)rsm0dr, (17) 

Jbot 

where Fg is the horizontal component of the local angular mo- 
mentum flux vector. Figure |9] shows the results from our solar 
model. The similarity between the the full model and the corre- 
sponding diagrams in Fig. 1 1 of Brun & Toomre (2002) is strik- 
ing, especially for their model AB, where the total horizontal 
transport is smallest. The L-effect is the most powerful trans- 
porter of angular momentum. It is balanced by the meridional 
flow and the viscous transport together, which are of compara- 
ble strength with the meridional flow being slightly stronger. 

Figure shows how the balance of the three transport 
mechanisms shifts as the rotation rate varies for the F8 star For 
fast rotation (4 d) the pattern is almost the same as for the solar- 
type star with the baroclinic term switched off. The meridional 
flow and the L-effect roughly cancel out each other. The vis- 
cous part of the stress tensor is of minor importance. The result 
for a period of 10 d looks very similar to the Fig.|9l The L-effect 



dommates and the How is stronger than the viscosity, though 
the latter is no longer negligible. For the solar rotation rate of 
27 d the viscous transport has become dominant. The L-effect 
and the meridional flow are of equal strength. Note that the con- 
tribution of the meridional flow has changed its sign. For very 
slow rotation (100 d) the horizontal L-effect is negligible. The 
meridional flow (driven by the radial shear) is balanced by the 
viscous transport (Kippenhahn 1963). 

5. Conclusions 

We have applied the mean-field theory of angular momentum 
transport to the convection zone of a main sequence star of 
spectral type F8. The only free parameter is the viscosity pa- 
rameter Ci, which is calibrated by applying the model to the 
case of solar differential rotation. For rotation periods in the 
observed range the resulting rotation pattern is similar to the 
solar one, though the horizontal shear is about twice as strong. 

The shear vanishes in both the limiting cases of very fast 
and very slow rotation. For fast rotation the Taylor-Proudman 
theorem requires the rotation rate as constant on cylindrical sur- 
faces. The only way to fulfill this constraint together with the 
boundary conditions is rigid rotation. For very slow rotation the 
horizontal L-effect vanishes and the rotation rate is a function 
of the radius only. 

The comparison of the contributions of the L-effect, tur- 
bulence viscosity, and meridional flow shows that a horizontal 
differential rotation can only be maintained if the L-effect is the 
dominant transporter To maintain a stationary rotation pattern, 
it must be balanced by the contributions from the flow and the 
viscous transport. Since it requires shear to be effective while 
the flow does not, it is also necessary that viscosity plays a ma- 
jor part in establishing the balance. A meridional flow driven 
by the shear alone will always counteract the force generating 
the shear. Thus, the angular momentum transport by the flow 
will be directed against that by the L-effect which generates the 
shear. For very large Taylor numbers the flow dominates over 
the viscous transport and the shear is strongly reduced. If, on 
the other hand, the flow is suppressed the angular momentum 
transport by the L-effect has to be balanced by the viscosity, 
which requires differential rotation to be effective. 

The solutions we find for the Sun and the F-star represent 
an intermediate state. The contributions from the flow and the 
viscous part of the Reynolds stress are of similar magnitude. 
Though not negligible, the meridional flow is not strong enough 
to balance the L-effect, leaving roughly half of the load to the 
viscous transport. Consequently, the model of Kiiker & Stix 
yields a value of 0.2 only for Sft/ft, rather than the 0.3 of 
the simpler model by Kiiker et al. (1993), where the Reynolds 
stress is the only transporter of angular momentum. We have 
therefore readjusted the model by choosing — 0.15 instead 
of 4/15, thus requiring more shear to maintain the same flux of 
angular momentum. 

The results suggest that the surface horizontal shear is 
mainly a function of the convection velocity, hence of the stel- 
lar luminosity. Our model does not produce a power law-type 
relation between rotation period and surface shear Local fits 
yield exponents that are greater than those found by observers. 
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Fig. 10. Horizontal flux of angular momentum for the F star with = 0.15 and rotation periods of 4 d (top left), 1 d (top right), 
27 d (bottom left), and 100 d (bottom right). The line styles are the same as in Fig.|9l 

but the scatter in the observations makes a direct comparison Strassmeier, K.G. 2002, Astron. Nachr., 323, 309 

impossible. Thompson, M.J., Christensen-Dalsgaard, J., Miesch, M., Toomre, J. 

2003, ARA&A, 41,599 
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